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Abstract. We show that Supersymmetric models with Type I seesaw neutrino masses support slow roll inflection point 
inflation. The inflaton is the D-flat direction labelled by the chiral invariant HLN composed of the Higgs(H), slepton(L) 
and conjugate sneutrino(N) superfields. The scale of inflation and fine tuning is set by the conjugate neutrino Majorana mass 
M v c ~ 10 — 10 12 GeV. The cubic term in the (quartic) inflaton potential is dominantly from superpotential (not soft Susy 
breaking) couplings. The tuning conditions are thus insensitive to soft supersymmetry breaking parameters and are generically 
much less stringent than for previous 'A-term' inflation scenarios controlled by mass scales ~ TeV . WMAP limits on the 
ratio of tensor to scalar perturbations limit the scale M controlling inflection point inflation: M < 7.9 x 10 13 GeV. 'Instant 

preheating' is operative and dumps the inflaton energy into MSSM modes giving a high reheat temperature : ~ Mj r 10 6 
GeV ~ 10 11 — 10 15 GeV. A large gravitino mass > 50 TeV is therefore required to avoid over closure by reheat produced 
gravitinos. 'Instant preheating' and NLH inflaton facilitate production of right handed neutrinos during inflaton decay and 
thus non-thermal leptogenesis in addition to thermal leptogenesis. We show that the embedding in the fully realistic New 
Minimal Supersymmetric SO(10) GUT requires use of the heaviest righthanded neutrino mass as the controlling scale but the 
possibility of a measurable tensor scalar perturbation ratio seems marginal. We examine the parametric difficulties remaining. 
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INTRODUCTION 

Primordial inflation is now generally accepted as the only viable mechanism for setting the initial conditions for 
Big bang Cosmogony in a way compatible with the extreme CMB homogeneity observed by increasingly accurate 
satellite maps [1] of the Microwave sky. The parameters of inflation (PR,n s , @k( n s)) measured so far can be accounted 
for by most of the many slow-roll inflation models proposed. With few exceptions these models use inflaton(s) that 
have no role to play outside of inflation since they have no connection with the known fields of Particle Physics. 
Such a connection is however necessary since the post-inflationary epoch must include reheating phases where the 
inflaton energy is converted into the matter and radiation observed today. Models driven by an inflaton composed 
of SM[2],MSSM [3] or GUT[4] fields thus carry an obvious appeal. Models of the second type type are typically 
based on slow roll inflation along "D-flat directions" in the MSSM field space and these are conveniently labelled by 
holomorphic gauge invariants formed from chiral superfields. Such models ( also called "A-Term Inflation" models 
[5, 6, 7]) typically require extreme fine tuning between the soft terms to ensure an inflection or saddle point of the field 
potential where the vacuum energy density drives a burst of inflation but nevertheless allows "graceful exit" due to the 
absence of a local minimum and the associated potential barrier which would prevent exit. Thus while they answer 
some of the relevant issues they have much scope for improvement. 

In [5, 6] an A-term inflation model was based on small neutrino yukawa couplings needed for realistic Dirac light 
neutrino masses. The inflaton field was a gauge invariant D-flat direction, NHL, where N is the right handed sneutrino, 
H is the MSSM Higgs doublet which gives masses to the up-type quarks, and L is the slepton field. When coupled 
with soft trilinear and bilinear supersymmetry breaking terms of mass scale ~ 100 GeV to 10 TeV the associated 
renormalizable inflaton potential can be fine tuned to achieve inflection point inflation consistent with WMAP 7 year 
data[5, 6]. 

The Type I seesaw[8]) mechanism offers a more attractive explanation for small neutrino(Majorana) masses 
(m v ~ (niy) 2 /M v c) based on large right handed neutrino masses M v c >> Ms- It is natural to ask if theories with 
supersymmetric Type I seesaw masses also support inflation. The popular Leptogenesis [9] scenario as well as the 
realistic Susy Minimal SO(10)GUT strongly hint at right handed neutrino masses in the range 10 6 to 10 12 GeV. So 
for Vb-l ~ Mx > 10 16 GeV the superpotential couplings fa,A= 1,2,3 which generate M v c ~ /aVb-l, are very small 
{Ja ~ 10~ 9 to 10~ 4 ) and can also give rise to a cubic term in the quartic inflaton potential. Thus the required ingredients 
for inflation are already present in Supersymmetric Type I seesaw models. 



Issues regarding natural values for superpotential couplings come into focus when viewed in the context of the 
so called Minimal Left Right supersymmetric models [10] and their embedding in GUT models[ll, 12]. SUSY Left- 
Right Models are advantaged due to their protection of R-parity as a gauged discrete symmetry, which provides a 
stable lightest supersymmetric particle (LSP) which has the properties required to serve as WIMP dark matter. They 
simultaneously and naturally implement Seesaw mechanisms for neutrino masses [10]. Moreover such models have 
also been incorporated in the realistic and predictive New Minimal Susy SO(10) grand unified theories(NMSGUT)[13, 
14] where all the hard parameters of the MSSM are fitted in terms of fundamental parameters of the GUT and 
soft SUSY breaking parameters (of the Non-Universal Higgs masses (NUHM) type) defined at the Unification scale 
Mx ~ 10 16 — 10 18 GeV. Such GUTs have viable Bino dark matter candidates and make distinctive predictions for the 
type of SUSY spectra observable at the LHC. In 2008, well before the discovery of Higgs mass of around 125 GeV in 
201 1-2012 and the consequent realization that a general framework such as the phenomenological MSSM (pMSSM) 
requires that the soft trilinear couplings A t j, be large, we concluded[13] that the NMSGUT would be falsified by its 
failure to fit the down type quark masses unless Ao,fi were in the tens/hundreds of TeV : leading to a mostly decoupled 
mini-split supersymmetry type superspectrum with only the LSP, gauginos and possibly a light slepton in the sub-TeV 
range. The experimental data has now forced this realization on practitioners of MSSM parametrology[15]. In the 
NMSGUT it was a prediction. 

Taken together with the possibility of small values for the light generation Yukawa Dirac couplings it is possible 
to implement viable inflection point inflation by suitable tuning at the supersymmetric level itself. This is technically 
more appealing than a tuning applied to soft susy parameters which, being unprotected by SUSY, are unstable. We 
examine the reheating dynamics briefly. We then derive derive the embedding of SSI in the NMSO(10)GUT and the 
necessary tuning conditions and show how to satisfy them explicitly. 



GENERIC PROPERTIES OF RENORMALIZABLE INFLECTION POINT INFLATION 

In this section we outline the essential features of inflection point inflation deriving from a quartic potential of a 
single complex scalar field (p. Since the angular degree of freedom has positive curvature and cannot support inflation 
one may assumed it fixed at its minimum by relaxation and focus on the remaining real field whose potential is 

genetically V = ^<p 4 - ^<j> 3 + ^-(j) 2 . 

The fine-tuning A = AM implies that one has a saddle point (V'(0 O ) = V"(0o) = 0) at the field value O = 
So one defines fine-tuning parameter A through A = AM^J 1 — A (A = j3 2 /4 in the notation of [6]). The inflection 
point(V"Oo) = 0) is at <fo = ^(1 - A + 0(A 2 )). For small A 

V ^ = y « = £( 1+4A ) : V\^) = a = ^± ; V"'(to) = y =^J(l-2A) (1) 

If h is tiny Vo >> M 4 and 0o >> M/h; / is small with h, while a is small by tuning. Large vacuum energy and flatness 
around 0o and starting near 0o with << 0q imply the universe executes slow roll inflation as rolls through an 
interval of width A0 <~ Vq/jM 2 below 0o- The standard slow roll parameters are defined as(M p = 2.43 x 10 18 GeV) 

M 2 n V" Ml Ml y' y n o Ml MtV'V" Mtay 

77(0) = ^^y(0-</>o) ; £ {<t>) = ^C- f^(a + l^-^f ?(^) ; | = ~ -^L (2) 

The observed CMB is a combined spectrum of modes which exited the horizon during inflation. We approximate 
it as a single spectrum from a representative("pivot") mode that exits the co-moving horizon when <j) = <Pcmb- This 
is the field value near <po where the inflation giving rise to observable effects today kicks in (when Ncmb e-folds of 
inflation are remaining). The power spectrum and spectral index we see today are then (Ncmb)) and « 4 (0(A^cmb)) 
respectively. The small first and third Taylor coefficients a,y determine [7, 16, 17] the measured parameters of 
inflation (PR,n s ) once the field values (<j>cMB,<l>end) at the time of horizon entry of the "pivot" momentum scale 
(kpivot — 0.002 Mpc~') and at termination of the slow roll are fixed[7, 17](on the basis of an overall cosmogonic 
scenario and the consistency of the slow roll approximation (77 (<j) en d) ~ 1) respectively). The observable number of e- 
folds Ncmb — N(<Pcmb) is the number of e-folds of inflation left to occur after <p crosses 0cmb (the field value when the 
representative primordial fluctuation length scale (l p i VOt — k p[ l m[ ) becomes larger than the comoving horizon ( 1 / a^HjJ). 
Plausible inflationary cosmogonies require 40 < Ncmb < 60 and this severely restricts the inflation exponents. 



The slow roll inflation formula for the power spectrum of the mode that is leaving the horizon when the inflaton 
rolls to and the corresponding spectral index and it's variation with momentum is([16]) 

W) = 247r 2^4 £W ;n,(0) = l+2Ti(0)-6e(0) ;^M = ^^ = -16 e 77+24 e 2 + 2^ (3) 

In practice e,£, are so small in the narrow region near 0o where slow-roll inflation occurs that their contribution to 
n s is negligible. 3>k{ n s) is negligible i.e. the spectral index is scale invariant in the observed range, as is allowed by 
observation so far. 

The field value at the end of inflation <j) eni i is defined by rj((j) mc i) ~ 1 which gives 0o — §end = —Sr. - In the slow roll 



approximation = —V'((j))/3H >> /H, where H = y V(0o) / (3A/|) is the (constant) inflation rate during slow roll 
inflation. One has the N — link (which can be exactly inverted[28], without assuming that <j) en( i << 0(iV) [7])): 

j'tyend H 2 Pi V i Ff _ Fy 
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A^pivo! is estimated using the standard Big Bang thermal cosmogony, giving [17] Af p i vo t = 65.5 + In rh M where 
p r h is the energy density after reheating and Vq the potential value during inflation. Due to rapid thermalization in 
this model(see below) the two are equal and then since the scale is set by V 1 / 4 ~ M/y/h one finds N p i VOt = 46 to 
55 = 5 1 ± 5 to be a reasonable estimate. 

To search for sets of potential parameters M,h, A compatible with the observed Pr,hs,Ncmb one uses the definitions 

Vo (n s -l) 
£cMB ~ lA^Mfa ; t]cmb 2 (5) 



and from these deduces OCcmb, §cmb using the eqns.(2) 

'()1]CMB 

yM 2 ; ~™»-v-™» Mp 2yM 4 



0cmb = 0o + ~rlr^ ; °fcMB = V^ £ cmbtj ^~ri£r (6) 



The required fine-tuning A is then 

llCCcMB , M 4 16/z 2 M p 2 

A= ^=W ( =^- (1 "^ )) (7) 

OcmBj A should emerge real and positive and using {(XcmbAcmb} in the formula for Ncmb one should obtain a sensible 
value in the range Ncmb = 5 1 ± 5. Using eqns.( 1 ,5) in eqn(4) we can solve accurately for the required relation between 
h,m,A using an interpolating function [28]. The result is that Ncmb ~ 50, Zq w solves the exact equations to a 
good approximation and one obtains the generic constraints : 

^3*^2.75x10-* M1 A ^14^ ^ 2 

M Mp^ MB ^cmb Ml N 2 CMB P R 

We then have viable inflation with inflaton energy and Hubble rate 

V ~ ^- - (M) 3 x 10 25 GeV~ 10 43 - 10 61 GeV 4 ; H - - 10 3 - !0 n GeV (9) 

The fine-tuning measure grows as M 2 so that j3 = can be as large as 10~ 2 forM <~ 10 12 GeV. In our scenario due 
to the large value of the inflaton mass parameter M ~ 10 6 to 10 13 GeV compared to M ~ TeV in the case of MSSM 
inflation[3] or Dirac neutrino inflation[5, 7] the fine-tuning required is quite mild and removes much of the motivation 
for complicated just so hybrid inflation scenarios. 

We can also estimate the ratio r of power in Tensor and Scalar CMB fluctuations using r = 2Vq/ (3k 2 PrMp). On 
using eqn.(8) 

M f (10) 



.95 x 10 13 GeV 

This makes the observation of tensor perturbations in such a scenario hard unless M is near its upper limit. 



SUPERSYMMETRIC SEESAW INFLATON MODEL 



In this section we introduce a toy one generation Supersymmetric seesaw inflation scenario model with gauge group 
SU(3)x St/ (2) x U(1)r x U (I)b-l that captures the essential features of our scenario. The essential fields beyond the 
MSSM are a right handed Neutrino chiral multiplet N[l ,1,-1/2,1] and a field S[l , 1 , 1 , —2] whose vev generates the 
large Majorana masses M v (10 6 — 10 14 GeV) for the conjugate neutrinos v£ = Na via a renormalizable superpotential 
coupling 3\/2fABSv^Vg. Additional superheavy fields £2, serve to fix the vev of S as in Minimal Supersymmetric Left 
Right Models (MSLRMs)[10] and in GUTs that embed them [11, 12, 13]. Neutrino Dirac Yukawa coupling is present 
in the superpotential : W = y v NLH + ... where L[l,2,0, — 1],//[1,2, 1/2,0] are the Lepton doublet and up type Higgs 
respectively. The relevant D-flat direction extends out of the minimum of the supersymmetric potential corresponding 
to the breaking of the gauge group down to the MSSM symmetry 

St/ (3) x SU (2) x U(l) R x U{\) B -l -> St/(3) x SU (2) x U{l) Y (11) 

This leads to a Type I seesaw plus MSSM (SIMSSM) effective theory. After symmetry breaking the MSSM hyper- 
charge Y — 2Tjr + (B — L) where T^r is the U(1)r generator. Unlike the case of the Dirac neutrino masses scenario 
[5] B — L is not a gauge symmetry down to low energies. This can have important consequences for nucleosynthesis 
and matter domination since the heavy right handed neutrinos must find a non-gauge channel to decay through. The 
flat-direction associated with the gauge invariant NLH is then specified as 

N = v = h = -^ = <j)e ie ; </>>0, 0e[O,27z;) (12) 
v3 

The additional fields £2, are coupled to S so that extremization of the SUSY potential using Fq = 0, D a \^=o = fixes 
the vev of S: < S >= o/^/l without constraining the inflaton field (p. This is as in the Minimal Susy LR models [10] 
and renormalizable Susy SO(10) GUTs [12, 13] which are our inspiration and target. 
At scales <j> ~ a » Ms where SUSY is exact the relevant superpotential is given by: 

W = 3 V3yNvh + 3/ V2SAW +...=y<p 3 + fVlSy 2 + ... (13) 

where h,f, a can be taken real without loss of generality. The equations of motion of the unperturbed vacuum imply 
< F$ >= 0, < S >= a j\[2. The right handed neutrino Majorana mass will be M v c = 6/(7. 
This superpotential leads to an inflaton potential 

V susy = \3y(p 2 + 2fd(p\ 2 + 2\f(p 2 \ 2 

= f 2 [(2 + 9y 2 )</> 4 +12y</) 3 <7cose+4a 2 </> 2 ] (14) 

Here y — y/f and fa sets the mass scale. Minimizing with respect to gives = % so we can focus on just the 
real part of <p and set <p = — <j> with <j> real and positive near the inflection point but free to fall into the well around 
<p = and oscillate around that value. In addition one also expects a contribution to the potential from the fi term for 
the Higgs doublets together with SUSY breaking quadratic and cubic soft terms, which we assume to be of the type 
generated by supergravity, but with non universal Higgs masses, i.e of the form: 

V soft = [A Q (y(p 3 +fV2S(p 2 )+h.c]+mj^\f\ 2 + m 2 H \H\ 2 + ml\H\ 2 

f 

= f 2 [yA ^acos3e+A Q a 2 ^ 2 cos2e + rfila 2 ^ 2 ] (15) 

here mo = mo/fa,Ao = 2Aq/ fa. The soft mass mo receives contributions from the sfermion and Higgs soft masses 
as well as the jj, term : m\ — (2m 2 ~ + mj 1 )/3. Here H are the sfermion and up type Higgs soft effective masses at 

the unification scale (m 2 ^ = m 2 H + |jll| 2 ). Since these masses and Aq should be in the range 10 2 - 10 5 GeV while the 
righthanded neutrino masses lie in the range 10 6 — 10 12 GeV, it is clear that mo,Ao are small parameters and even 
for the large values of mo,Ao ~ 10 s GeV found in the NMSGUT mo,Ao << 1. Thus these terms cannot significantly 
change 9 = n assumed earlier. The total inflaton potential is then 



V tot = f 2 ((2 + 9y 2 )<j> 4 - (A + U)ya^ + (A + m 2 +4)a V) . 



(16) 



Thus we have a generic quartic inflaton potential of the same type as in Section 2 but the parameter values in the case 
of Type I seesaw are quite different from the light Dirac neutrino case. We have the identification of parameters 



A = 



fyj 12(2 + 9f) 
3f(A +12)ya 



v/(2 + 9y 2 ) 
M 2 = 2f 2 a 2 (4+A + ml) 



A 2 

(1--^—) 
1 16M 2> 



9f(A + \2) 2 , 



32(2 + 9y 2 )(A +m 2 +4) 

For seesaw models the natural magnitude for the neutrino Dirac mass is, > lMeV (i.e \y%\ > 10~ 5 and then 
the limit m v « O.OleV for the lightest neutrino (assuming direct hierarchy) implies M v c > 10 6 GeV). Since the 
preferred values for the Susy breaking scale are smaller than 100 TeV (at most) it follows that the maximum value 
of |Ao|, |mo| ~ 0.1 and they could be much smaller for more typical larger values of the conjugate neutrino masses 
M v c ~ 10 8 to 10 12 GeV. It is then clear from the corresponding range A ~ 10~ 12 to 10~ 4 that the coupling ratio y — y/f 
becomes ever closer to exactly y = 4/3 as M increases and even forM <~ 10 6 GeV differs from 1.333 only at the second 
decimal place. Thus to a good approximation h — 6y/6f. Then it follows that 

/ ~ lO-^iO.lV^) . M „ 10 -25.38 ± 0.17 ( ^_ )2 . A „ 10 -78.93 ± 0.47 ( ^_ )4 (lg) 

The range M ~ 10 6 6 to 10 10 6 GeV corresponds nicely to 10 16 GeV < a < 10 18 GeV: as is natural in single scale Susy 
SO(10) GUTsfll, 12, 13, 14]. / increases with a with values below 10~ u achievable in the NMSGUT only with 
difficulty. In MSLRMs, since there are no GUT constraints on <7, one may assume somewhat wider ranges for these 
parameters. 

In all relevant cases A < 10~ 4 is required. Thus the above equations imply that y 2 must be close to the value 

2 _ 64 4+Ap + m 2 
y ° 9 16-8A -32m 2 +A 2 k ' 

Here Aq , mo ~ O (Ms / M v c ) < < 1 , hence yo is rather close to 4/3 and the equality is very close for larger M <~ fa since 
then Ao,mo are tiny. The measure of severity of fine tuning j3 = VA ~ 10 -2 — 10~ 6 compares quite favourably with 
the case of the MSSM or Dirac neutrino inflaton since there j3 <~ 10 -12 to 10~ 10 due to the low values of the inflaton 
mass in those cases. The dominant component of the fine tuning in the present case is a fine-tuning of superpotential 
parameters, which is radiatively stable due to non renormalization theorems. Specially for large a > 10 16 GeV the 
Type I Susy seesaw can provide a rather attractive inflationary seesaw with a natural explanation for neutrino masses 
and weaker tuning demands on the radiatively unstable Susy breaking parameters than the extreme and unstable fine- 
tunings demanded by typical inflection point scenarios and in particular the Dirac neutrino model [5] . Moreover, unlike 
the chaotic sneutrino inflaton scenario [18, 19], no trans-Planckian vevs are invoked. 



Reheating and Leptogenesis 

The post inflationary dynamics of our model bears an intimate relation to previous studies of models with 'instant 
preheating' mechnism[20] and specially the MSSM flat direction inflection point inflation model[21] and preheating 
model[22] with strong coupling to the MSSM Higgs. Supersymmetric seesaw inflation offers an attractive synthesis 
precisely fulfilling the need expressed in [22] : 

"There have been many models of leptogenesis. A hallmark of our model is the economy of fields. The only 
undiscovered fields are the inflaton, <j), the standard model Higgs, h, and the right-handed neutrino, N. There are 
very good reasons for suspecting that all exist! The only unfamiliar aspect of our model is the strong coupling of the 



inflaton field to the Higgs field. While there is no reason to preclude such a coupling, it would be very interesting to 
find particle-physics models with a motivation for the coupling. " 

Due to the gauge(H,L) and third generation yukawa(H) coupled components of the inflaton the inflaton energy will 
decay very rapidly (with decay time Tdec « Hr 1 ^ ~ (hM p )/M 2 ) by the 'instant preheating" mechanism[20, 22, 21]. 

Thus the reheating temperature r r /, ~ T max ~ V^ 4 ~ M/h l l 2 <~ 10 11 — 10 15 GeV. The parametric dependence is 
identical to that found in [21]. The difference in scales arises only because the inflaton mass M~ 10 6 -10 12 GeV in our 
model is much larger than the inflaton mass parameter m§ <~ 0.1 — 10 TeV in [21] coming from soft Supersymmetry 
breaking. 

In preheating ("% type") degrees of freedom, with masses(m z <~ gty{t)) and decay rates (r ~ g 3 (j)(t)) proportional 
to <j>(t), are produced non-perturbatively every time the inflaton field crosses zero. This happens because the % modes 
are ultra-light for a sufficiently large time interval around the zero crossing time during which adiabaticity is violated ( 

> (0% : where (Ok is the oscillation frequency at wave number k). Here % modes are the components of the H,L, u c l ,ul 
chiral superfields and the W±,B gauge superfields. They can be identified as the fields which become massive given 
background values of the three components of the inflaton (V, v[,h°). Then with the usual superpotential 

W = y u Q L Hu L L + y d Q L Hd L L + y v LHN + y l LHe L L + ... (20) 

y leads to massive ml,m£; y v leads to massive e/Xone combination of the three e^ ), h°,h + ,Vl, v[; y l leads to massive 
h~,e c L (one combination). Since < H,N,L> preserve U(l) em , the gauge couplings give masses to Z (which forms 
a Dirac supermultiplet with (v —h~o)/\/2) and W± (form a pair of Dirac supermultiplets with l~,h+). The inflaton 
vev leaves the down quark and gluon/gluino fields and Jiq, and some combinations of the l£,l^ fields with light 
(MSSM type) masses. These light (y/-type) fields will form the first step in the decays of the % field. As < > 
again increases after crossing zero the % modes become heavy and unstable and as a result decay rapidly(within a 
time x dec ~ j^-j « m^ 1 ) to the light (mostly coloured) MSSM modes to which they are coupled strongly coupled. 

A fraction ~ 10 _1 of the inflaton condensate energy passes into the light MSSM modes with every crossing resulting 
in complete transfer within <~ 10 2 oscillation times. x osc ~ 1 << H^j ln ~ (hM p )T osc /M ~ (1 — 150)t ojc . Once 
the energy is in the light(y/) modes MSSM interactions rapidly complete thermalization. Rapid decay of the inflaton 
oscillation amplitude leaves the light modes to thermalize the energy dumped by the inflaton into a radiation bath of 
all modes: which are no longer ever heavy because the inflaton has decayed. The reheating temperature is 

T rh ~ (^-) 17 X 1/4 - T nlax ~ 10 11 - !0 15 GeV (21) 

where g* = 228.75 is the effective number of MSSM degrees of freedom. This reheating temperature is well above 
that required to produce relativistic populations of gravitinos : which are unacceptable if their lifetimes are larger than 
the nucleosynthesis time Tjv ~ 1 sec since their decay after nucleosynthesis would destroy the created nucleons. The 
straightforward and generic resolution of this gravitino problem is if the graviton masses are sufficiently large so that 
the gravitinos decay before nucleosynthesis[23] : T grav ~ 10 5 sec(^^-) 3 « %n ~ 1 sec. Thus Supersymmetric seesaw 

Inflation also indicates that the scale of supersymmetry breaking -as indicated by the gravitino mass- should be above 
50 TeV; as is also found by fitting of fermion data in the NMSGUT[ 13] . Large reheat temperatures also ensure abundant 
thermal production of righthanded neutrinos after inflation. Their CP violating decays into leptons can drive thermal 
lepto-genesis [9] for generating the observed baryon to entropy density hb/s ~ 10~ 10 . Non-thermal leptogenesis is 
also possible [22] since the Higgs field H is itself a % type field and coupled to the righthanded neutrinos. During 
inflaton oscillations the Higgs mass m/, <~ g2§ fluctuates below and above M v c ~ fa. CP violating Higgs-righthanded 
Neutrino inter-conversion [22] leads to (non-thermal) Leptogenesis which will add to the thermal leptogenesis. The 
complication in the present case that the L,H and N components of the inflaton have different decay rates implies a 
proper analysis must track the separate evolution of all three fields making up the inflaton using the equation of motion 
and Boltzmann equation for the relevant degrees of freedom. This requires a separate numerical study to expose the 
interplay of the couplings /a, yAB, 82- The study of this evolution and the operation of Leptogenesis in these models is 
now in progress. 

INFLATION AND NEUTRINO MASSES IN THE NMSGUT 

Finally we embed SSI in the New Minimal SO(10) GUT (NMSO(IO)GUT or NMSGUT). The NMSGUT is a realistic 
Susy SO(10) model[13, 14, 25] that successfully fits the known fermion mass-mixing data in terms of GUT parameters 



and provides structural reasons for suppression of the dangerous operator dimension d = 4,5 Baryon violation typical 
in Susy GUTs[14]. It furthermore makes distinctive predictions of a mini-split supersymmetry spectrum made viable 
by large A,/i-terms and with a characteristic normal s-hierarchy. Neutrino flavour plays a key role in enabling 
NMSGUT inflation : the inflaton is composed of third generation conjugate sneutrino, first generation left slepton 
(sneutrino) and the 7\r — 1/2 Higgs. 

The NMSGUT Higgs field vevs {210(a),/?, a), 126(a)} = ii,l26(S = a) break 50(10) Gmssm while preserving 
Supersymmetry at Mx- An explicit Susy preserving solution of symmetry breaking in terms of a cubic equation for a 
complex variable x and depending on a single parameter ratio ^ was found by us[12]. The mass spectra implied[12, 13] 
by this analytic solution for the the MSGUT vacuum are the basis of our detailed Renormalization Group and threshold 
effect analysis[12, 13, 14]. Inclusion of threshold corrections raises the unification scale close to the Planck scale and 
can lower the gauge coupling at unification. We shall use the notation and results of [26, 27, 13, 14]. 

To embed SSI corresponding to a NLH type flat direction we show there is a corresponding flat direction of the full 
GUT potential which rolls out of the MSGUT minimum (that has the SIMSSM as its effective theory). The relevant 
fields are the GUT scale vev fields Q. = {(0,p,a,w, g},S — G and the (6) possible components hj, hf,i = 1...6 of the light 
MSSM Higgs doublet pair H,// together with the chiral lepton fields L A ,v c A ,A = 1,2,3. The relevant superpotential is 
then[12, 13] 

W = 2V2(h AB h x - 2V3 f AB h 2 - 8 ab(Ii5 + iV3h 6 )) + h T Jf(< ft >)h + 4V2f AB av A v B + Wh(ft) (22) 

where 

W n (ft, a) = m(p 2 + 3a 2 + 6co 2 ) + 21 (a 3 + 3p(0 2 ) + (M+ri(p + 3a- 6(o))o6 (23) 

and 

-^^■|fc,v,L=0 = -^\h,v,L=0 = D a{&)\h,v,L=0 = (24) 

here h AB , g AB , f AB are the yukawa coupling matrices of the three matter 16-plets to the 10,120,126 Higgs multiplets 
respectively. Jf? is the Higgs doublet mass matrix[27, 12, 13]. Equation (24) defines the MSGUT vacuum[12]. 

Of the 5 diagonal D-terms of SO(10) only those corresponding to the generators T^l, 
Tj,r,B — L are relevant for vevs ft, a and out of equilibrium inflaton mode composed of v,V c ,Iiq. The vevs 
ft, a do not contribute to these D terms or cancel so 

Oil = f (-I>o| 2 + E|v A | 2 ) 

L i=l A 

D3R = f (£|M 2 -2|/mo| 2 -E|v a | 2 ) 

L i=\ A 

Db-l = \Jlgufc(\vA\-\vA\ 2 )+2\h 4Q \ 2 ) (25) 

where only h 4a = has B - L = +2, T 3R = - 1/2 and thus Y = 1 while all others have T 3R = 1 /2 and B - L = 0. 
The D-flatness conditions are 

E|va| 2 = E|/* 2 |-E|v a | 2 + 2|/mo| 2 (26) 

A i A 

In MSGUTs the MSSM Higgs doublet pair is defined by fine tuning Det(J^) ~ so that its lightest eigenvalue 
/I — M w ~ 1 TeV specifies the jU term in the superpotential of the SIMSSM : W = fiHH + .... The doublet pair H,H 
is a linear combination [12, 27, 26] of the 6 doublet pairs of the the NMSGUT : 

hi^UijHj hi = UijHj (27) 

where U, U diagonalize the doublet mass matrix Jff : U T JifU — Diag{/J. , ....,Af¥} to positive masses. They are 
calculated with jU = = Det(Jf). The so called Higgs fractions : a, = Ua,0Cj = Ua , determine the grand unified 
formulae[12, 13] for the SIMSSM fermion yukawas. For tree level yukawa couplings replace -4- OiH, a,H. For 
example the neutrino Dirac coupling is {{h AB ,g AB J AB )=2\ f 2(h AB ,g AB J AB ) ) 



y v AB = h AB a { - 2V3 f AB a 2 - g AB {a 5 + iV3a 6 ) 



(28) 



From the V = \Fj r \ 2 only the light Higgs doublet H can contribute. To get small yukawas the involvement of the 
lightest generation is unavoidable. Thus we take Va = Vi. Taking Va = Vj the tuning constraint is at best of form 
|yn| 2 <~ 10(|y2i| 2 + |^3i | 2 ) : this is impossible to satisfy with normal neutrino yukawa coupling hierarchy. Choosing 
Va = v| as the conjugate neutrino component of the inflaton is more helpful in satisfying the fine tuning condition. 
Thus our inflaton ansatz is 

*'=75 h = ^f^M> (29) 

Note how the Higgs fraction 0:4 enters as T = 1 — 2 1 0:4 1 2 . It happens the solutions we have found earlier [13] can have 
I OC4 j ■ — ' 0.5. It is not inconceivable that T ~ is achievable without destroying the realistic fermion fits to the fermion 
data. 

We use generic Supergravity(SUGRY)-NUHM generated soft terms in terms of a common trilinear parameter Aq 
but different soft mass parameters m 2 pih 2 hj for the 16 plets and the different Higgs. Repeating the analysis of section 3 
with the NMSGUT superpotential and the new ansatz we obtain the parameter identifications 

h = ^[(/V)n+r(|/i3i| 2 +4||3i| 2 + (>'V t )33)+4|/33| 2 r 2 )] i 

A = i(16|/33||y3 V ill^|Vf + 4b3 V i|VfA o cos(30 ( j-20,v i )) 
M 2 = — |/33| 2 |a| 2 r+-Ao/ 3 3|a|rcos(30 (f -20 > ,v)+2mg (30) 
The fine tuning condition A = AM is now 

b3Vl|2= 9A rf -8A„"(l+r) [l>-nl 2 + l^| 2 + r(fei| 2 +4|g3i| 2 + b 3 V 2l 2 + b3 V 3l 2 )+4|/33| 2 r 2 ] (31) 



Where 



K = 1 + ^2-008(30,-2^)-' 3 ™ 2 ° 



4M 3 y \6M 2 T 
A d = (1 + ^-^(30^-20^ )) 2 (32) 

and M3 = /33 1 a\ M v >>M$ imples A n j are both very close to unity. So as before the fine tuning condition is essentially 
between hard parameters as in GUTs and in sharp contrast to MSSM inflaton models[3]: 

In NMSGUT fits the strong hierarchy ^33 1 >> \y^\ » \yn \ » \yn \ > \yu\ holds . So one must tune 

TwO i.e \a 4 \~-^= (34) 
v2 

to a good accuracy. The MSSM doublet H is almost exactly 50% derived from the doublet in the 210 plet ! The yukawa 
tuning condition is only 

l^i I 2 = mi\ 2 +\yli\ 2 ) (35) 

which is easy to enforce in the NMSGUT. 

There is an additional demand coming fromeqn(8) : h 2 /Mi ~ (y v ^y v )ii/Mj, ~ 10~ 25 which is, at first glance, much 
harder to enforce. However [14] large wave function corrections [29] to the GUT(Fj-"' £ ')-MSSM(}y) yukawa coupling 
relation due to the circulation of heavy fields within loops on the lines entering the yukawa vertex imply : 

Y f = (1 + Aj) • {Yf)tree ■ (1 + Af) ( 1 + A H ± ) (36) 



TABLE 1. Illustrative example of relevant parameters from an accurate fit of the fermion spectrum 
in the NMSGUT which is compatible with inflationary scenario. All masses are in GeV. Xx,z are the 
accuracies of the fits to 18 known fermion mass/mixing parameters at Mx t z- 



Parameter 


Value 


Parameter 


Value 


Xx 


0.4458 


M h o 


122.99 


Xz 


0.1426 


Mx 


7.08 x 10 17 


h 


1.066 x 10" 3 


f\Jl 


2.59 x 10~ 8 , 4.405 x 10~ 5 


h 


2.44 x 10~ 4 


A„ 


0.999999 


M 


3.043 x 10 11 


A^ 


0.999999 


r 


4.343 x 10~ 5 


^tuning 


0.989 


1*1 


4.69 x 10 15 


M X 


5.25 x 10 17 


A Q (M x ),m Q (Mx) 


-5.235 x 10 5 , 1.260 x 10 4 


H,B(Mx) 


4.316 x 10 5 , -1.128 x 10 11 


M 2 - 

\*H \!\Aa \ 
Mf 


-1.498 x 10 11 


M H 


-1.448 x 10 11 


50.254,63.930 


\Ot\ 


0.707 


4.86 x 10 13 




1.181 x 10 9 ,2.01 x 10 12 


\yl\tree\ 


1.997 x 10~ 4 


\y2\tree\'\y\\tree\ 


4.489 x 10~ 5 , 1.640 x 10~ 6 


Log w {h 2 /M) 


-18.706 


Volenti 


3.579 x 10 52 , 2.153 x 10 15 


N pivot, NqmB 


54.22,4.78 x 10" 4 


A,J3 


8.82 x 10~ 12 ,5.92x lO" 6 



Due to the large number of heavy fields the dressing of the Higgs fields can be rather large (<~ 10 2 ). We earlier 
calculated [14] the dressing for the 10-plet component of the MSSM Higgs. Above we showed that a completely 
independent line of argument requires that the doublet H be 50% derived from the 210-plet. Thus the lengthy 
calculation of the wave function corrections for each of the six GUT doublets contributing to the MSSM doublet 
is necessary. Even from the partial calculation [14] we see that the large value of the wave function dressing makes 
the GUT tree level matter fermion yukawa couplings (i.e {fiAB,gAB,fAB}tree an d therefore all the (y^ B ) tre e ) required 
to match the SIMSSM couplings at much smaller than they would be without these corrections. This has the 
important consequence of suppressing d = 5 B-violation operators since they depend on these yukawas and have no 
Higgs line. Since it is the tree level couplings that enter the formulae for the infiaton dynamics in the full GUT it is 
easier to satisfy eqn.(8). Because of this and the relatively large value of M <~ M3 it should be be possible to achieve 
the required fine tuning using the full wave function dressing [24]. 

Embedding in the GUT has overturned our naive assumption that the lowest intermediate scale would govern 
inflation. Instead it is rather the largest. While setting us the problem of finding solutions to the tuning condition, 
compatible both with an accurate fit of fermion masses and acceptable values of inflationary power spectrum and 
spectral index, it emphatically shows that the soft terms have little role to play in the fine tuning which belongs rather 
to the GUT and intermediate scale physics only. Thus the physics of SIMSSM driven inflation is in sharp contrast to 
the Dirac neutrino mass driven inflation[5, 6] in the MSSM extended by U{\)b-l and right handed neutrinos. Our 
analysis makes it clear that they lie counterpoised not only as regards the nature of neutrino mass but also as regards 
the nature of inflation and its regulating mass scale besides their degree of naturalness. 

An example of the relevant parameters from an accurate fit of the complete fermion spectrum in the NMSGUT which 
has also been tuned to make it as compatible as possible with the inflationary scenario presented is seen in Table 1 . 
More details may be found in [28]. The fine tuning between the yukawas proceeds as anticipated with 1 — Y = 1 = A n c i. 
The remaining problem is that h 2 /M ~ 10~ 19 GeV is too large. As a result the number of e-folds Ncmb is much smaller 
than required. However as explained the formulae used seriously underestimate the Higgs wave function corrections. 
Search of the huge parameter space has just begun. We may well hope for a completely realistic fit compatible with 
inflation in due course. The detailed analysis of reheating behaviour and Lepto-genesis in this model is also underway. 
Since the fits of fermion masses also yield a value for the CP violation parameter relevant for Lepto-genesis the 
incorporation of the SSI scenario in the NMSGUT may eventually yield useful additional constraints which will serve 
to narrow the parameter space further. 
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